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CLAIMS 



1-27. (Cancelled) 

28. (New) A method of determining a position estimate based on an updated Kalman filter, 
comprising: 



receiving a first measurement Li based on a first signal with wavelength Xi and 
frequency fi; 

receiving a second measurement L 2 based on a second signal with wavelength X 2 
and frequency f 2 ; 

selecting a model a of distance dependent and distance independent errors in the 
first and second measurements, wherein the model a is selected from a = 
Xi/X 2 , a = X 2 /X u and a = 1 ; 

based on the model a, calculating a double differenced variance matrix: 



D(L X ,L 2 ) = 
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K+K ef W ref ... W ref 

Kef w 2 +w ref ... w ref 




D 22 =(R 2 L2 +a 2 R 2 Ll -2aR UL2 ) : ! \ ! 

Kef Kef - K-2+Kef 

Kef Kef - Kef 




wherein 



wherein 



R L 2=(^[RL-^l2+(RL +RL p +RD-B 2 ])"; 



wherein 



Rn,L2 =tt( r L +RL P +RL)-(^kRL +Rl P +RD-B 2 ; 



wherein Rnon-dist.n is a non-distance dependant measurement error for the 
Li measurement; 

wherein R„ 0 n-dist,L2 is a non-distance dependant measurement error for the 

L 2 measurement; 
wherein Ri on is measurement error due to ionospheric delay; 
wherein Rtrop is measurement error due to tropospheric delay; 
wherein Rqa is measurement error due to orbit bias; 
wherein each W is: W n = 1 .0 + 7.5e^ /15 ; 

wherein n is a satellite ordinal of a plurality of satellites; 
wherein W re f corresponds to a reference satellite of the plurality of 
satellites; 

wherein E is the elevation angle of each satellite of the plurality of 
satellites; 

updating a Kalman filter with the calculated variance matrix; and 
determining a present position estimate based on the updated Kalman filter. 
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29. (New) The method of claim 28, wherein the updated Kalman filter with the 
calculated variance matrix is K = D{L l ,L 2 ) new Hl{H T 2 D(L l ,L 2 ) new H T 2 + V 2 Y l ; 

wherein H 2 is a transpose of a design matrix H 2 for the measurement L 2 ; 

wherein V 2 is a residual measurement of measurement L 2 ; 

wherein D{L x ,L 2 ) nm! = D{L l ,L 2 ) old - K old H 2 D(L X , L 2 ) old ; 

wherein D{L^L 2 ) nexv is the calculated variance matrix including the received 
measurements Li and L 2 ; 

wherein D(L x ,L 2 ) old is a previous calculated variance matrix which was 

calculated before receiving both of the received measurements Li and L 2 ; 
and 

wherein K old is a previous Kalman filter which was calculated before receiving 
both of the received measurements Li and L 2 . 

(New) A computer readable medium having computer-executable instructions, 

receiving a first measurement LI based on a first signal with wavelength W and 
frequency fl; 

receiving a second measurement L2 based on a second signal with wavelength X2 
and frequency f2; 



30. 

comprising: 
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selecting a model a of distance dependent and distance independent errors in the 
first and second measurements, wherein the model a is selected from a = 
MA2,a = A2A,l,anda=l; 

based on the model a, calculating a double differenced variance matrix: 



D{L X ,L 2 )- 



I Ai 0 

o a, 
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wherein 



R L2 ={-^{RL- dist ,L2 + (RL +Rl P +Rl b )-B 2 ]y; 



wherein 
Rn 



1 



" ( R L + R L P + Rorb ) • ("4 Rln + Rlop + R orb ) 



wherein R n0 n-dist,Li is a non-distance dependant measurement error for the 
Li measurement; 
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wherein R n on-dist,L2 is a non-distance dependant measurement error for the 

L 2 measurement; 
wherein Ri on is measurement error due to ionospheric delay; 
wherein Rtrop is measurement error due to tropospheric delay; 
wherein Ro rb is measurement error due to orbit bias; 
wherein each W is: W n = 1 .0 + 1.5e~ Ells ; 

wherein n is a satellite ordinal of a plurality of satellites; 
wherein W re f corresponds to a reference satellite of the plurality of 
satellites; 

wherein E is the elevation angle of each satellite of the plurality of 
satellites; 

updating a Kalman filter with the calculated variance matrix; and 
determining a present position estimate based on the updated Kalman filter. 

3 1 . (New) The computer readable medium of claim 3 1 , wherein the updated Kalman 
filter with the calculated variance matrix is: 

K = D(L l ,L 2 ) new H T 2 {H T 2 D(L l ,L 2 ) new H T 2 +V 2 y l ; 

wherein Hi is a transpose of a design matrix H 2 for the measurement L 2 ; 

wherein V 2 is a residual measurement of measurement L 2 ; 

wherein D(L y ,L 2 ) new = D(L„L 2 ) old - K oU H 2 D{L l , L 2 ) old ; 

wherein D(Z 1 ,£ 2 ) new is the calculated variance matrix including the received 
measurements Li and L 2 ; 
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wherein D{L x ,L 7 ) old is a previous calculated variance matrix which was 

calculated before receiving both of the received measurements Li and L 2 ; 
and 

wherein K oM is a previous Kalman filter which was calculated before receiving 
both of the received measurements Li and L 2 . 
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